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ABSTRACT 

It is known that merging galaxies with luminous starbursts and high 
far-infrared luminosities tend to have higher i?i_o = 12 CO(J=1-0)/ 13 CO(J=1-0) 
integrated line intensity ratios (i?i_o — 20 — 50) than normal spiral galaxies 
(i?i_o — 5 — 15). Comparing far-infrared luminosities [L(FIR)] with those of 
12 CO(J=1-0) and 13 CO(J=1-0) for a sample of normal and starburst galaxies, 
Taniguchi & Ohyama found that the observed high i?i_o values for the luminous 
starburst mergers are attributed to their lower 13 CO line intensities by a factor 
of 3 on the average. They suggested the following two possibilities; in the 
luminous starburst mergers, 1) 13 CO is underabundant with respect to 12 CO, or 
2) exitation and/or optical depth effects are responsible for the change in Ri-o- 
In this paper, we investigate the second possibility using higher transition data 
of both 12 CO and 13 CO emission lines. Applying the same method proposed 
by Taniguchi & Ohyama to both 12 CO(J=2-l) and 13 CO(J=2-l), we find that 
13 CO(J=2-l) is also depressed with respect to 12 CO(J=2-l). This suggests 
that the 13 CO gas may be underabundant in the high-i?i_o starburst mergers 
although we cannot rule out the possibility that excitation and optical depth 
effects are still affecting i? 2 -i, for example due to the large velocity widths in 
the CO emission lines. Additional observations of both 12 CO and 13 CO lines at 
J > 3 are required to better constrain the conditions of the molecular gas in 
luminous starburst galaxies. 

Subject headings: galaxies: emission lines - galaxies: starburst - interstellar: 
molecules 
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1. INTRODUCTION 

Since starburst phenomena affect both the chemical evolution of galaxies and the 
physical conditions of the interstellar medium in galaxies, an understanding of starburst 
activity is one of the important issues in astrophysics. Starbursts occur in dense molecular 
gas clouds which are generally located in the central regions of galaxies. Massive stars 
formed in the starbursts provide a negative feedback to the parent and ambient molecular 
gas clouds in the form of intense radiation fields and strong stellar winds as well as 
subsequent supernova explosions. Therefore it is important to investigate the molecular gas 
properties of starburst galaxies (Young & Scoville 1991; Henkel, Mauersberger, & Baan 
1991; Sanders & Mirabel 1996 and references therein). 

One of the interesting properties of molecular gas in starburst galaxies is that galaxy 
mergers with luminous starbursts (i.e., L(FIR) ^ 1O 11 L : hereafter luminous starburst 
mergers) tend to have higher i?!_ (= J[ 12 CO( J=1-0)]/7[ 13 CO( J=l-0)]) integrated line 
intensity ratios than normal spiral galaxies (Aalto et al. 1991, 1995, 1997; Casoli et al. 
1991; Casoli, Dupraz, & Combes 1992a, 1992b; Hurt & Turner 1991; Turner & Hurt 1992; 
Garay, Mardones, & Mirabel 1993; Henkel & Mauersberger 1993; Henkel et al. 1998). One 
possible explanation for the higher R±-o values is that they are due to the inflow of disk gas 
with high 12 C/ 13 C abundance ratios, possibly combined with a 12 C enhancement caused by 
nucleosynthesis in massive stars (e.g., Henkel et al. 1998). However, recently, Taniguchi 
& Ohyama (1998a; hereafter T098) have compared far-infrared luminosities [L(FIR)] 
with CO luminosities, L[ 12 CO( J =1-0)] and L[ 13 CO( J =1-0)], for a sample of normal and 
starburst galaxies, and have found that the observed higher i?i_o values are associated 
almost exclusively in the luminous starburst mergers and appear to be attributed to a 
lower intensity of 13 CO(J=1-0) with respect to 12 CO(J=1-0) when compared to normal 
galaxies. T098 suggested either that 13 CO is underabundant with respect to 12 CO, or that 
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the 13 CO(1-0) level population is more depressed relative to 12 CO(1-0) due to excitation 
and/or optical depth effects, leading to the high i?i_o in the luminous starburst mergers 
studied in their paper. In this paper, we investigate the second possibility using available 
higher transition data of both 12 CO and 13 CO emission lines. 

2. 12 CO(J=2-l)/ 13 CO(J=2-l) INTEGRATED LINE INTENSITY RATIO 

T098 demonstrated that the comparison of L(FIR) with both L[ 12 CO] and L[ 13 CO] 
provides a powerful tool for understanding the origin of the high-i?i_ in luminous starburst 
mergers (see also Taniguchi & Ohyama 1998b). If the observed high i?i_o values are 
attributed to exication and optical depth effects, for example higher gas kinetic temperatures 
and/or denser gas clouds, then this could possibly be discerned in the measured values of 
_! or, if necessary, in even higher transition line ratios. In order to examine if this is the 
case, we investigate the excitation properties of both 12 CO and 13 CO molecules for a sample 
of normal and starburst galaxies using data currently available from the literature. 

We have compiled 12 CO(J=2-l) and 13 CO(J=2-l) integrated intensities from Aalto 
et al. (1995) and Casoli et al. (1992b). The integrated intensity ratio J[ 12 CO(J=2- 
l)]//[ 13 CO(J=2-l)] is referred as Ri-\- Our sample consists of 24 galaxies and includes 
objects with extreme infrared luminosities such as the ultraluminous infrared galaxy Arp 
220. These integreated intensities are then used to compute CO luminosities; L(CO) is 
defined as L(CO) = A x /(CO) K km s _1 pc 2 where A is the observed area in units of pc 2 
and 7(CO)= jTljri~ l dv K km s _1 where T£ is the observed antenna temperature corrected 
for atmospheric extinction and rj is the main beam efficiency. 

The FIR data are compiled from the IRAS Faint Source Catalog (Moshir et al. 1992). 
The FIR luminosities are estimated using L(FIR)= 4tt j D 2 1.26 x Kr n [2.58 x S{60) + 5(100)] 
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(ergs s -1 ) where 5(60) and ,5(100) are the IRAS 60 /im and 100 /im fluxes in units of Jy and 
D is the distance (Helou, Soifer, & Rowan- Robinson 1985). Distances of nearby galaxies 
are taken from the Nearby Galaxies Catalog (Tully 1988); distances of other galaxies are 
estimated using a Hubble constant H = 75 km s' 1 Mpc -1 with Vgsr (recession velocity 
with respect to the Galactic Standard of Rest) given in de Vaucouleurs et al. (1991). The 
compiled data are given in Table 1. All of the data presented here have been corected 
for beam-size (see Aalto et al. 1995; Casoli et al. 1992b). Although our sample is not 
statistically complete, it is the largest sample compiled so far. T098 defined the class of 
high-i?i_ galaxies by adopting the criterion of Ri-o > 20. Using this limit, the present 
sample contains the following seven high-i?^ objects; NGC 1614, NGC 3256, NGC 4194, 
NGC 6240, Arp 220, Arp 299, and IRAS 18293-3413. 

In Figure 1, we compare L[ 12 CO( J=2-l)] with L[ 13 CO(J=2-l)]. For reference, we 
also show the comparison between L[ 12 CO( J=l-0)] and L[ 13 CO( J=l-0)] in the left panel 
which is taken from T098. Although two high-i?i_ galaxies, NGC 4194 and NGC 6240, 
have significantly lower L[ 13 CO( J=2-l)] with respect to L[ 12 CO( J=2-l)], the remaining 
galaxies have R2-1 ratios within the upper range (i.e. ~10-30) found for i?i_o- 

In Figure 2, we compare L(FIR) with both L[ 12 CO( J=2-l)] and L[ 13 CO( J=2-l)]. 
L[ 12 CO( J=2-l)] appears to be correlated with L(FIR) (in an integrated intensity versus 
flux plot, i.e. after removing the the D 2 effect from Figure 2, there is a good correlation), 
but the scatter is larger than observed for the correlation between L[ 12 CO( J=l-0)] and 
L(FIR) (T098). On the other hand, the correlation between L(FIR) and L[ 13 CO( J =2-1)] 
is poorer than that between L(FIR) and L[ 12 CO( J=2-l)] because the majority of the 
high-i?i_ galaxies have lower L[ 13 CO( J=2-l)] by a factor of 3 than what would be expected 
from the correlation for the normal-i?i_o galaxies. Thus, we find that both 12 CO(J=2-l) 
and 13 CO(J=2-l) show similar behavior as observed in the J=l-0 transition (T098). In 
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order to show this more clearly, we present a diagram of Ri-o versus R2-1 for our sample 
in Figure 3. In Table 2, we summarize the statistical properties for both the high-i?!_ and 
the "normal- i?i-o" galaxies. 

3. 12 CO(J=3-2)/ 13 CO(J=3-2) INTEGRATED LINE INTENSITY RATIO 

Published measurements of extragalactic 13 CO(J=3-2) emission are available only for 
M82 (Tilanus et al. 1991; Wild et al. 1992) and IC 342 (Wall & Jaffe 1990). Tilanus et al. 
(1991) obtained 13 CO(J=3-2) spectra at three positions in M82; the center and the two 
peaks of the circumnuclear ring located ±12 arcsec from the center using the 15 m JCMT 
with a beam size of 14 arcsec (FWHM). They obtained R3-2 — 15 for the nucleus while 
~ 10 for the circumnuclear ring. On the other hand, Wild et al. (1992) measured both 
12 CO(J=3-2) and 13 CO(J=3-2) emission lines using the IRAM 30 m radio telescope in 
February 1992 and obtained both J[ 12 CO( J=3-2)] = 1334 K km s' 1 at (Act, AS) = (-5, 
-5) and J[ 13 CO( J=3-2)] = 70.4 K km s" 1 at (Aa, AS) = (-7, -5) where Aa and A5 are 
offsets from the nucleus position in right ascension and declination, respectively, in units 
of arcsec. Although the measured positions are slightly different, these measurements give 
an integrated intensity ratio, R3-2 — 18.9. This together with the results by Tilanus et al. 
(1991) indicates that R3-2 in the nuclear region is ~ 15 - 20. Since this value is nearly 
the same as the threshold value which defines the class of high-i?i_o galaxies (T098), it 
suggests that the J = 3 transition is still not high enough to allow dissentanglement of 
radiative transfer and abundance effects in M82. 

As suggested by Taniguchi & Ohyama (1998b), it is possible that a large value of R can 
be attributed to the effect of superwind activity (i.e., the possible destruction of dense gas 
as well as dust grains, and the large velocity widths observed in the CO outflow). In fact, 
in M82, R3-2 is higher in the nuclear region than in the starburst ring (Tilanus et al. 1991). 
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Since M82 is indeed a superwind-starburst galaxy (Bland & Tully 1988), the higher -R 3 _ 2 
value in the nuclear region of M82 could possible be due to superwind activity. However, 
Tilanus et al. (1991) suggested that intensity ratios of the three lowest transition lines of 
12 CO and 13 CO can be explained if 13 CO is overabundant with respect to 12 CO just like 
what has been measured in the Galactic center, being contrary to our interpretation. 

Another measurement of the 13 CO(J=3-2) line was obtained at the central region of 
IC 342; R 3 ^ 2 ^ 7.7 (Wall & Jaffe 1990). However, since the beam size of 13 CO(J=3-2) 
observation (24 arcsec) is different from that of 12 CO(J=3-2) one (15 arcsec), the above 
value may not be reliable. 

Until more extragalactic CO(J=3 - 2) data is obtained it is clearly impossible to draw 
any firm conclusions about R3-2 in luminous starburst mergers as well as normal galaxies. 

4. DISCUSSION 

We have shown that the R 2 -i ratio is also high (typically > 20) in the high-_R!_ 
galaxies. Furthermore, the R3-2 value in the nuclear region of M82 suggests that this ratio 
may also be high in luminous starbursts, but this is only for one object. However, it is 
noted that the 12 CO(J=3-2)/ 12 CO(J=l-0) integrated intensity ratio of starburst galaxies 
is often found to be higher than that in normal galaxies (e.g., Devereux et al. 1994 and 
references therein). Further, some nearby starburst galaxies such as M82 are detected in 
CO(J=4-3) (Giisten et al. 1993) and in CO(J=6-5) (Harris et al. 1992). The detection 
of these higher-transition CO lines suggests the presence of warm and dense gas clouds in 
starburst galaxies. 

Although the kinetic gas temperature is not necessarily comparable to the dust 
temperature, it is interesting to compare molecular gas properties with the dust 
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temperature which is measured from the IRAS 60 /im to 100 /im flux ratio, 5(60) /S( 100). 
In Figure 4, we compare i? 2 _i with S(Q0)/S(100) for the galaxies studied here. We also 
compare Ri-o with S'(60)/S'(100) for reference (the left panel). Though no tight correlation 
can be seen in either of these diagrams, we find that both i?i_o and -R2-1 tend to increase 
with increasing S'(60)/S'(100). This tendency suggests that the galaxies with higher dust 
temperatures have higher R values on the average. In the lower panels of Figure 4, 
we show comparisons of S(60)/S(100) with both L[ 12 CO( J=2-l)]/L[ 12 CO( J =1-0)] and 
L[ 13 CO(J=2-1)]/L[ 13 CO(J=1-0)] ratios. Interestingly, we find no correlation in both the 
diagrams. Average ratios of both L[ 12 CO( J=2-l)]/L[ 12 CO( J=l-0)] and L[ 13 CO(J=2- 
1)]/L[ 13 CO(J=1-0)] are 0.81 ±0.17 and 1.15 ±0.71, respectively for the high R^ starburst 
mergers. Therefore, there seems to be no significant difference in the excitation toward J=2 
between 12 CO and 13 CO. 

We also investigate whether a correlation exists between i?i_o and the luminosity 
ratio L[ 12 CO(J=2-1)]/L[ 12 CO(J=1-0)] for our sample (Figure 5). Since the high-i?i_ 
mergers tend to have higher dust temperatures (see Table 2), it is likely that their CO 
kinetic temperatures are also higher than those of the normal- i?i_ galaxies. These high 
temperatures could lead to both higher i?i_ and to higher L[ 12 CO( J=2-l)]/L[ 12 CO( J=l- 
0)] as demonstrated by large velocity gradient models (Sakamoto et al. 1994, 1997). 
However, there is no such tendency as shown in Figure 5. 

Finally, we investigate whether there is any relationship between i?i_o, R2-1 and the 
L(FIR)/L(CO) ratio which is generally considered to give a measure of star formation 
efficiency. Figure 6 shows that the high-i?i_ mergers as a group have a higher mean 
L(FIR)/L(CO) ratio as well as a higher mean R2-1 ratio than the "normal-i?i_ " galaxies. 
But, other than that there appears to be no clear correlation between either i?i_o or 
R2-1 with respect to L(FIR)/L(CO). Therefore, we are led once again to suggest that 
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"superwinds" may be the best explanation for what produces the high-i?!_ values. Indeed, 
it should be noted that nearly all of the high-i?!_ galaxies show morphological and/or 
spectroscopic evidence for superwinds (Taniguchi & Ohyama 1998b). Both the abnormally 
large velocity gradients assocated with these superwinds and the possible destruction of 
dense gas clouds by the dynamical effect of the superwind activity, could possibly combine 
to reduce the observed intensity of the much more optically thin 13 CO in the lower- J 
transitions relative to the 12 CO emission. Further tests of this hypothesis will require 
measurements of the J > 3 transitions of both 12 CO and 13 CO for the galaxies in Table 1. 

We would like to thank Seiichi Sakamoto for useful discussions. YO was supported by 
a Grant-in- Aid for JSPS Fellows by the Ministry of Education, Science, Sports and Culture. 
This work was supported in part by the Ministry of Education, Science, Sports and Culture 
in Japan under Grant Nos. 07055044, 10044052, and 10304013. 
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Fig. 1.— Diagram of L[ 12 CO(J=2-l)] vs. L[ 13 CO(J=2-l)] (right panel). For reference, we 
also show a diagram of L[ 12 CO( J=l-0)] vs. L[ 13 CO( J=l-0)] taken from T098 (left panel). 
The high-i? galaxies (i?i_o > 20) are shown by filled circles. 

Fig. 2.— Diagrams of R 2 -i (top), L[ 12 CO(J=2-l)] (middle), and L[ 13 CO(J=2-l)] (bottom) 
against L(FIR). The symbols have the same meaning as those in Figure 1. Alphabets in the 
upper panel mean; a: NGC 6240, b: NGC 4194, c: Arp 299, d: NGC 3256, e: Arp 220, and 
f: IRAS 18293-3413. 

Fig. 3. — Diagram of i? 2 — i vs. Ri-o- The symbols have the same meaning as those in Figure 
1. 

Fig. 4.— Diagrams of R 2 -i (upper right), i?i_ (upper left), 12 CO(J=2-l) to 12 CO(J=1-0) 
luminosity ratio (lower right), and 13 CO(J=2-l) to 13 CO(J=1-0) luminosity ratio (lower 
left) against S (60) / S (100) . Note that R 2 -i = 140 for NGC 6240. The data of i?i_ are 
taken from T098. Alphabets in the upper left panel mean; a: NGC 4194, b: NGC 6240, c: 
NGC 3256, d: NGC 3256, e: NGC 1614, f: NGC 3256, g: IRAS 18293-3413, h: NGC 3256, 
i: Arp 220, j: Arp 299, k: Arp 299, 1: IRAS 18293-3413, m: ESO 541-IG 23, and n: Arp 
220. 

Fig. 5. — Diagram of R ± _ vs. 12 CO(J=2-l) to 12 CO(J=1-0) integrated intensity ratio. 
The symbols have the same meaning as those in Figure 1. The large velocity gradient models 
taken from Sakamoto et al. (1997) are show for the case of X (CO) / dv / dr = 10~ 5 [km s _1 
pc -1 ] -1 where X(CO) is the fractional abundance of 12 CO and dv/dr is the velocity gradient. 
Kinetic temperatures and densities of molecular gas are labeled. 

Fig. 6. — Diagrams of Ri-o and i? 2 -i against L(FIR)/L(CO) ratio. The symbols have the 
same meaning as those in Figure 1. Note that the data shown in the left panel are taken 
from T098. 
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Table 1. Properties of the sample of galaxies 
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O 1 


o.Zl 


l .16 


I .61 




U.OOO 


lU.oo 


in i 

lz.1 




l.U 




1 


N6240 


98 a 


9.72 


8.08 


9.32 


7.18 


0.816 


11.52 


43.6 


140 


0.8 


0.3 


2 


N6951 


24 


8.40 




8.36 




0.353 


10.21 






0.9 




1 


N7552 


20 


8.72 


7.56 


8.70 


7.74 


0.710 


10.63 


14.2 


9.1 


1.2 


1.9 


1 


N7582 


18 


8.50 


7.28 


8.58 




0.673 


10.39 


16.7 




1.1 




1 


IC2554 


17 


8.09 


6.97 


7.48 




0.475 


9.95 


13.0 




1.1 




1 


U2855 


20 


8.52 


7.44 


8.44 


7.46 


0.387 


10.40 


11.9 


9.5 


1.1 


1.4 


1 


U2866 


16 c 






7.32 


6.19 


0.792 


10.19 




13.5 






1 


Circinus 


4 d 


7.85 


6.83 


7.60 


6.40 


0.787° 


9.78° 


10.5 


15.9 


1.1 


1.2 


1 


Arp 220 


74 a 


9.72 


< 8.39 


9.42 


8.16 






> 21.8 


18.2 


0.7 


> 1.0 


2 


Arp 299 


42 a 


8.99 


7.68 


8.50 


7.04 






20.7 


28.9 


0.9 


0.7 


2 


I18293 f 


73 c 


9.79 


8.56 


9.66 


8.70 


0.688 


11.46 


16.8 


9.1 


1.8 


1.4 


1 



Note. - (1) Name, (2) Distance (Mpc) taken from Tully (1988), (3) log L[ 12 CO(J = 1-0)] (K km 
s- 1 pc 2 ), (4) log L[ 13 CO(J =1-0)] (K km s" 1 pc 2 ), (5) log i[ 12 CO(J = 2-1)] (K km s" 1 pc 2 ), 
(6) log L[ 13 CO(J = 2-1)] (K km s" 1 pc 2 ), (7) S(60)/S(100), (8) L(FIR) (L ), (9) i?i_ , (10) R 2 -i, 
(11) i[ 12 CO(J = 2 - l)]/L[ 12 CO(J = 1 - 0)], (12) L[ 13 CO(J = 2 - l)]/L[ 13 CO(J = 1 - 0)], and (13) 
Reference. 

a Estimated using a Hubble constant H = 75 km s _1 Mpc -1 with Vgsr given in de Vaucouleurs et 
al. (1991). b Sanders et al. (1995) c Strauss et al. (1992) d Freeman et al. (1977) c Lonsdale et al. (1989) 
f I18293 = IRAS 18293-3413 



References. — (1) Aalto et al. 1995; (2) Casoli et al. 1992 
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Table 2. Comparisons of physical properties between the high-i?!_ mergers and the normali? 1 _ galaxies 





Units 


High-i?i_ mergers 


Normal-i?i_o galaxies 


Ri-o 




31.5±12.7 


12.2±3.2 


R2-1 




47.3±45.2 


12.7±4.9 


log L[ 12 CO (J =1-0)] 


(Le) 


9.30±0.39 


8.50±0.63 


log L[ 13 CO(J=l-0)j 


(£0) 


7.87±0.34' 


7.44±0.52 


log L(FIR) 


(L Q ) 


11.36±0.34 


10.26±0.53 


log L(FIR)/L[ 12 CO(J=1-0)] 




2.06±0.27 


1.76±0.18 


T d 


(K) 


44.4±2.6 


38.2±3.6 



